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Abstract. We calculate the Hilbert action for the Bondi-Sachs metrics. It yields the 
^ | Einstein vacuum equations in a closed form. Following the Dirac approach to constrained 

O systems we investigate the related Hamiltonian formulation. 
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cn ; 1 Introduction 

This paper is a folow up of Refs [T| [2] in which we studied the Bondi-Sachs metrics [3] jlj 
in relation to the Penrose conformal approach to the asymptotical flatness [H| jB] • Here we 
consider the Lagrangian and Hamiltonian formulation of the Einstein theory (in vacuum) 

CT) • for these metrics. We derive the exact Einstein equations and we single out those of them 

which are constraints from the point of view of an evolution with respect to the Bondi 
time u. We perform also a preliminary study of the related canonical formulation in a 

"q \ manner similar to that of Goldberg [2j [H] (see also p] , where the Ashtekar variables are 

used for a null foliation). Our approach differs substantially from canonical formalism 
on asymptotically null hyperboloids (see e.g. [Ej [H] and references therein). 

Unlike in PQ[2], m this paper we consider Bondi-Sachs metrics using the affine dis- 

.£h ! tance instead of the luminosity one. Thus, we assume that spacetime metric g of signa- 
ture H takes the form 

g = du(g 00 du + 2dr + 2ujAdx A ) + r 2 gABdx A dx B . (1) 

The coordinates u and r are interpreted, respectively, as a retarted time and a distance 
from matter sources and x A , A = 2, 3, are coordinates on the 2-dimensional sphere 
S2 with the standard metric sab- The function g 00 , the form u)Adx A and the tensor 
gAsdx A dx B are regular on [a, b] x [r , 00] x S 2 - In order to guarantee the asymptotical 
flatness of g at null infinity these objects are assumed to admit the following expansions 
with respect to Q — r _1 

ffoo = 1 - 2Mfi + 0(n 2 ) , (2) 
Co A = ijA + ttn A + 0(tt 2 ) , (3) 
9ab = s ab + Vtn AB + 0(tt 2 ) , (4) 
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with coefficients depending on u and x A . In accordance with the Einstein equations we 
also assume that 

s AB n AB = . (5) 

Above properties of g assure that the conformally equivalent metric g = Q 2 g is well 
defined up to the scri J? + [Hj, where Q = 0. 

Equation © is a reminiscent of the luminosity condition 

det g AB = det s AB (6) 

considered in original papers of Bondi et al [3] and Sachs In this paper we replace 
©by © and 

fib = 1 (7) 

in order to simplify the Lagrangian formulation. Now r is the affine parameter for null 
geodesies generating surfaces u =const, 

g aP u\ a d p = 8 r , (8) 

but M is not the Bondi mass aspect. Under assumptions (JIJ)-© one can always pass 
from (0) to © by means of the transformation 

r = ar' , (9) 

where 

a = (^i)V4 = i + bQ 2 + 0(fi 3 ) (10) 
det g 2 

and b does not depend on r. Under this transformation the fields n AB and ip A do not 
change and 

M' = M-b„, (11) 



k' a = k a + b iA . (12) 

Note that M' is the Bondi mass aspect. 

Substituting ©-© into the Einstein vacuum equations yields 

IpA = ~\ nB A\B , ( 13 ) 

9ab = -(1 + \^ 2 n 2 )s AB + nn AB + 0(fi 3 ) (14) 
o 

(hence b = j^n 2 ) and the following relations between the u-derivatives, denoted by a 
dot, of the coefficients M, n AB and k a 

& = -\n AB n AB , (15) 

o 
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*>A = {\nh - ^M) >A - ^n AB h BC ]c + {-h AC n c B - ^n c [AlB]c ) lB . (16) 

Here 

M = M- l -nh + -n AB ^ AB , (17) 

n 2 denotes n A Bn AB and the metric — s A b is used to lower and raise the indices A, B 
and to define the covariant derivatives \a (note that this convention changes in the next 
sections). 

Concerning a physical meaning of different fields on we note that M can be used 
instead of the Bondi mass aspect to describe the total energy related to a section E of 
the scri [12] [2j 

E(E) = — I Mda . (18) 

Equation (JT5j) is responsible for the energy loss formula [Hj |3] ■ The traceless tensor 
fiAB corresponds to the Bondi news function. Functions ka play an important role in a 
description of the angular momentum at null infinity (see e.g. [H] and references therein). 



2 Lagrangian and equations 



In this section we consider the Hilbert action (modulo the constant c 3 /16ttG) 

I = - J d 4 x^/\g~\R (19) 

for metrics Since the metric components are fixed the variational principle for 
yields all the Einstein equations except G lM = 0, where G^ v is the Einstein tensor 
of g. The lacking equations are also fulfilled provided all other Einstein equations are 
satisfied as well as asymptotic conditions ©-© and equations (|I5jL (fTB^l (see jl] [H] 
and discussion at the end of this section). We postpone here the problem how to obtain 
the latter equations from a Lagrangian defined on J^ + . We treat these equations as 
definitions of functions M and ka in terms of their initial values and functions uab- 

In order to express the action ()19)) in a form corresponding to the ADM action ^H] 
we procceed as follows. First we complete the metric by the term —edr 2 , where e is 
a positive parameter. Then we find the ADM action related to the spacelike foliation 
by u =const and, after removing singular terms, we take its limit when e — > 0. This 
approach leads to the following representation of the action integral I (up to boundary 
terms) 

I = Jd 4 X d[t AB PAB - KP\a + PABP AB ) - 2u A pJ B - ^AB^J] . (20) 

Here x l = r, a = \J det <jab aud all operations on indices A, B, including covariant 
derivatives, are defined by c/ab- Function A is related to g 00 via 

g 00 = \ + £o A u A (21) 
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and functions p AB are given by 

Pab = -\gAB,i + (}na) tl gAB ■ (22) 

In order to control behaviour of fields at null infinity it is convenient to express the 
action (J2"U|) in terms of the unphysical metric g = Q 2 g 

g = du[(X + u A u A )du - 2dVt + 2uj A dx A ] + g AB dx A dx B . (23) 

Then one obtains 

I = J d 4 x a[g AB p AB + \{p% - p AB p AB ) - 2u A pJ B - \g A B^ B ] , (24) 
where now x 1 = Q = -, all operations on indices A, B correspond to the metric g A B, 

Pab = ^gAB,i - Qna),ig AB (25) 

and cr = y/det g A s- We notice that u A = u A , A = Q 2 \, g AB = Q 2 ()ab an d p AB = 
p AB — VLg AB . Conditions © and © are equivalent to 

X = Q 2 -2MQ 3 + 0(fi 4 ), (26) 



u A = n 2 ^ A + n 3 K A + o(n 4 ). (27) 

The asymptotic behaviour of g AB is governed by (jU) and (JSJ). 

Let us consider the action J as a function of A, uj a and g AB . The variable A appears 
in as a Lagrange multiplier, in a similar way as the lapse function in the ADM 
formalism. The variation of I with respect to A yields the equation 

P A aa -PabP AB = (28) 

which is equivalent to the Einstein equation Gu = 0. Equation (}2"%|) is a constraint from 
the point of view of the evolution with respect to u . In terms of g^ u it reads 

\gAB,ig A ^-\^<r),ii = ^ (29) 

Equation (J2~9"j) reduces the number of independent degrees of freedom in g AB to two. 
The variation of (|24|) with respect to uj a leads to constraints 

{ogABW^ih ~ 2vpJ B = (30) 

which are equivalent to the Einstein equations G\ A = (we keep a instead of aQ~ 2 
in order to shorten notation). It follows from (jHUj) that fields u A can be defined by a 
double integration of some functions of g A s and their derivatives tangent to the surfaces 
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of constant u. In this procedure the coefficients ka appear as the integration constants 
(hence, condition (fTBj) can be imposed). 

The variation of ()24)1 with respect to g AB yields dynamical equations equivalent to 
Gab = 0. They read 

9ab,i + ^9ab + ^9ab,i + 9 CD 9c{a,i(^b)\d - 9b)d ~ ^9b)d,i) ~ \{^ C 9ab,i)\c 

+U^9ab,i - pu{a\b)),i - \g A c9BDU C iUJ D i + [%S + \(g C D - 1u c \d)9 CD \ 
-\ gcD ^ x uj D x + u c (\na), lc + U™ C \c),i ~ i(A^),n - Ri - HUqab = 0, (31) 

where S denotes the l.h.s. of 

Under conditions ©-© expanding equations (j3~Uj) and (}3*T]) into powers of fl implies 
relations (pH?|) and (fT^jl . It follows from them and from equations (|TK|) and (fT7)|) that 
Q~ 2 G 0fl = on y + (equivalently, Vt~ 2 G 1 ^ = on J? + ). By virtue of the Bianchi iden- 
tities these boundary conditions assure that G lfl = everywhere. Thus, the variational 
principle for (fTTJj) yields a complete set of the Einstein equations under assumptions (fTB]) 
and (Jill). 

3 Hamiltonian 

Let us consider a canonical formulation based on the Lagrangian density 

C = d[g AB PAB + \{p A A>1 - PABP AB ) - 2u A pJ B - \g A Btu jcuf] (32) 

corresponding to the action ()24|). Since A is a Lagrange multiplier we take the space 
of fields (g AB , uj a ) as the configuration space. Elements of the phase space are given 
by (g AB , u A , ttab, IIa), where tiab and Ua are the momenta conjugate, respectively, to 
g AB and u A . Since the Lagrangian density C is linear in g AB the inverse Legendre 
transformation is not well defined. Following the standard procedure of Dirac JH] 
we treat equations 

KAB ~ OPAB = (33) 

as constraints which should follow from a total Hamiltonian. Apart from that, since L 
contains no time derivative of u A we must add further constraints 

n A = 0. (34) 

Constraints (jHHJ), (|3"4*|) together with (f2Hj) are primary constraints of the theory. Taking 
them into account leads to the following completion of the canonical hamiltonian density 
ir A B9 AB -C 

H = a[2u A p AB lB + ^g AB ujjuj B } - XaS + [i AB <\>ab + ^ A IT A , (35) 

where fi AB , v A are Lagrange multipliers and 4>ab denotes the l.h.s. of constraint (JHHJ) . We 
notice that Ti is not a sum of constraints, unlike in the ADM calculus. The Hamiltonian 
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functional is denned as H = J cPxH (we disregard here boundary terms and problems 
with possible internal singularities of surfaces u =const). 

The next step of Dirac's algorithm is to assure a proper propagation of all constraints. 
From 

{<P AB ,H} = Q (36) 

we get complicated equations which do not involve time derivatives and can be treated as 
conditions imposed on the multipliers \i AB and A. By virtue of the Hamilton equations 

9^ = ^ = ^, " A = ^ = " A (37) 
dir AB dll A 

equations (JSfiJl become equivalent to (}3*Tj) (using this fact one can easily reconstruct an 
explicit form of (jSEJ))- Thus, the dynamical equations (pTTj) are encoded in the conditions 
on the Lagrange multipliers. 

The propagation of the constraint ()29)1 gives another condition on \i AB 

H AB i9ab,i - Vab,x9 Ab a + 2//% = > ( 38 ) 

whereas constraint leads to equation (fHU|) . Since (f3T)|) does not contain any Lagrange 
multiplier it is a secondary constraint of the theory. Its time derivative yields equations 



\\°9abv B - au B (n AB + \^9ab)],\ + {y-AB,x ~ Hb9ac,i + ^9ab,i) IB 



\9bc,i^ BC \ A ~ H,iA - fi(\na) M = 0, (39) 



(here /i = fj, A ) which can be considered as conditions on the multipliers v A . 

Thus, we end up with conditions (j3Hj). (J3~%j) . (JHH1) on the Lagrange multipliers and 
with constraints (}29j) . (}30j) . (}33j) . (}34j) on the fields g AB , to A and their momenta tc ab , ^a- 
All these constraints are second class in Dirac's terminology (the lack of first class 
constraints is related to fixing coordinates in metric (^)). The nonvanishing Poisson 
brackets of the constraints are given by 

{4>ab, 4>(a CD )} = aa AB ,i + \olab°a + \aag ABA - oa c ^ A g B )c,\ 

-\\ooF D g C D,\ + cr« i + \ocas\gAB (40) 

{4>ab, S(f3)} = -\9ab(3,ii + \P,\9ab : i + \P{9ab,xi - g CD gAc,i9BD,i) (41) 
{^ A ,4>(a BC )} = -\[vu B ( a AB + \ocg AB )\,i + (a AB ,i - a B g AC> \ + \ag A B,\i B 

-\9BC,ia BC \ A ~ a M ~ a(\na) M (42) 

{U A ,Hl B )} = -^9ABl B 1 ),i. (43) 



Here ip A denotes the l.h.s. of (}3T)j) and a AB , [3, •y A are smooth smearing functions (e.g. 
?/>(7 A ) = J d 3 xipAl A and a = ot A A ). 

An alternative possibility of the Hamiltonian formulation is to treat the variables 
uj a in the Lagrangian (J32|) as functions of gAB defined by equations (}30|) (the resulting 
expressions for u A are still simpler than those for v A in the previous formulation) . The 
corresponding Euler-Lagrange equations remain equivalent to the Einstein equations. 
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Now the phase space consists of elements (g AB ,ttab)- The total Hamiltonian density 
is given by (jnU) with v A = 0. Equations flU and flM)) are the primary constraints. 
Their propagation yields equations (}36|) and (}38|) which define the Lagrange multipliers 
A and fi AB . In this formulation the number of constraints and Lagrange multipliers is 
substantially reduced but still the Hamiltonian picture remains quite complicated. 

4 Concluding remarks 

We transformed the Hilbert action for the Bondi-Sachs metrics to the forms ()20j) and 
()24|) similar to the ADM action. The variational principle for the action (J2U yields 
equations (j2T?j) -(f3~T ]) which are equivalent to the Einstein vacuum equations provided the 
r-independent equations (fi3j) . (fTBj) are satisfied. Equations (HHJ), (JHUj) are constraints 
from the point of view of the evolution with respect to the Bondi time u. They can be 
solved in terms of integrals over Q (or r). 

We analyzed briefly the corresponding Hamiltonian formulation. The Dirac method 
for constrained systems leads to the Hamiltonian (J3l)j) (or its truncation with v A = 0) 
with the Lagrange multipliers which carry the main information about dynamics of the 
system. A striking result is that the Hamiltonian does not vanish when the Einstein 
equations are satisfied. The Lagrange multipliers are defined implicitly via a complicated 
system of differential equations. There are six second class constraints and no first class 
constraints. 
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